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1.  INTRODUCTION 


Accurate  dynamic  models  are  required  to  establish  the  dynamic  response  of 
complex  satellite  structures.  Unfortunately,  analytical  dynamic  models  of 
complex  structures  agree  closely  with  properly  measured  mode  data  only  in  the 
first  few  modes.  The  effect  of  this  deficiency  can  be  minimized  by  using,  as 
the  dynamic  model,  the  measured  modes  directly.  This  approach  has  been  used 
successfully  on  numerous  satellite  programs.  Another  approach,  generally 
referred  to  as  system  identification,  is  to  adjust  the  analytical  dynamic 
model  in  an  attempt  to  improve  correlation  between  analytical  and  empirical 
modes.  A  relatively  large  quantity  of  work  has  been  published  in  this  field. 
However,  no  procedure  developed  to  date  has  achieved  wide  acceptance  in  the 
community. 

Numerous  goals  and  approaches  to  the  identification  problem  are  presented 
in  the  literature.  For  example^  Rodden^  published  a  procedure  for  estab¬ 
lishing  structural  inf luence' coefficients  from  modes  of  an  effectively  uncon- 

2  .  , 

strained  structure.  Subsequently,  Hall,  using  optimization  theory,  formu¬ 
lated  a  procedure  that  established  a  stiffness  matrix  such  that  the  resulting 
analytical  modes  matched  selected  empirical  modes  in  a  least  squares  sense. 
The  mass  matrix  was  assumed  to  be  exact.  A  procedure  to  adjust  an  analytical 

mass  matrix,  and  derive  an  approximate  stiffness  matrix,  was  published  by 

3  4 

Berman  and  Flannelly.  Ross  introduced  a  procedure  for  deriving  both 

^Rodden,  W.P.,  "A  Method  for  Deriving  Structural  Influence  Coefficients 
from  Ground  Vibration  Tests,"  AIAA  J. ,  Vol.  5,  No.  5,  May  1967. 

^Hall,  B.M. ,  "Linear  Estimation  of  Structural  Parameters  from  Dynamic  Test 
Data,"  AIAA/ASME  11th  Structures,  Structural  Dynamics,  and  Materials 
Conference,  Denver,  Col.,  22-24  Apr.  1979. 

^Berman,  A.  and  W.G.  Flannelly,  "Theory  of  Incomplete  Models  of  Dynamic 
Structures,"  AIAA  J. ,  Vol.  9,  No.  8,  Aug.  1971. 

■‘Ross,  R.C.,  Jr.,  "Synthesis  of  Stiffness  and  Mass  Matrices  from  Experimental 
Vibration  Modes,"  SAE  paper  710787,  Sept.  1971. 
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Che  mass  and  stiffness  matrices  from  measured  natural  frequencies  and  a  square 
modal  matrix  composed  of  measured  modes  supplemented  by  arbitrary  linearly 
independent  vectors.  A  similar  concept  was  recently  published  by  Zak,"’  who 
supplements  the  measured  mode  data  with  information  from  analytically  predicted 
modes . 


'Iterative  procedures  that  employ  statistical  parameter  estimation  to 
adjust  analytical  models  have  been  published  by  Collins,  et  al.^  and  Lee  and 
Hasselman.^  These  procedures  will  preserve  the  physical  configuration  of 
the  model.  However,  as  demonstrated  in  Ref.  7,  the  agreement  between  the 
adjusted  model  and  test  data  can,  for  some  parameters,  be  worse  than  the 
agreement  prior  to  adjustment. 


Identification  procedures  based  on  matrix  perturbation  theory  have  been 

8  9  10 

proposed  by  Chen  and  Mada,  Chen  and  Garba,  and  Chen,  et  al.  These  pro¬ 
cedures  assume  Chat  the  analytical  model  parameters  are  close  to  Che  correct 
quantities.  This  is  a  relatively  severe  restriction  since,  for  complex  struc¬ 
tures,  experience  indicates  that  at  least  a  few  parameters  will  hot  satisfy 
this  requirement. 


5zak,  M. ,  "Discrete  Model  Improvement  by  Eigenvector  Updating,"  ASCE  J. 
Engineering  Mechanics,  Vol  109,  No.  6,  Dec.  1983. 

^Collins,  J.D.,  et  al.,  "Statistical  Identification  of  Structures,"  AIAA 
J. ,  Vol.  12,  No.  2,  Feb.  1974. 

^Lee,  L.T.  and  T.K.  Hasselman,  "Dynamic  Model  Verification  of  Large 
Structural  Systems,"  SAE  paper  781047,  Nov,  1978. 

^Chen,  J.C.  and  B.K.  Wada,  "Criteria  for  Analyses-Test  Correlation  of 
Structural  Dynamic  Systems,"  J.  of  Applied  Mechanics,  June  1975. 

^Chen,  J.C.  and  J.A.  Garba,  "Analytical  Model  Improvement  Using  Modal  Test 
Results,"  AIAA  J. ,  Vol.  18,  No.  6,  June  1980. 

^*^Chen,  J.C.,  et  al.,  "Direct  Structural  Parameter  Identification  bv  Modal 
Test  Results,"  AIAA  paper  83-0812,  AIAA/ASME/ASCE/AHS  24^*^  Structures, 
Structural  Dynamics  and  Materials  Conference,  Lake  Tahoe,  Nev. ,  2-4  May  19S3. 
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Identification  procedures  have  also  been  developed  using  constrained 

minimization  theory.  Baruch  and  Itzhack  introduced  formulations  to  adjust 

11  12 

analytical  stiffness  and  flexibility  matrices  such  that  the  resulting 

analytical  dynamic  model  modes  are  identical  to  the  analytically  orthog- 

onalized  test  modes  used  in  the  identification.  These  procedures  assume  that 

13  . 

the  mass  matrix  is  correct.  Berman  introduced  a  formulation  that  modifies 

the  mass  matrix  and  assumes  that  the  measured  modes  are  exact.  Subsequently, 
14 

Berman  and  Nagy  combined  Che  mass  matrix  adjustment  procedure  of  Ref.  13 
with  the  stiffness  matrix  adjustment  procedure  of  Ref.  11  to  establish  Che 
so-called  analytical  model  improvement  (AMI)  procedure.  This  procedure  will 
yield  a  model  whose  modes  agree  exactly  with  those  used  in  Che  identifica¬ 
tion.  However,  as  demonstrated  in  Ref.  10,  Che  analytical  mass  and  stiffness 
matrices  can  be  dramatically  altered.  Particularly  troublesome  is  the  modi¬ 
fication  of  stiffness  coefficients  from  values  of  zero  to  large  magnitude 
nonzero  values.  Clearly,  the  introduction  of  load  paths  that  do  not  exist  in 
Che  actual  hardware  is  undesirable. 

As  can  be  ascertained  from  Che  preceding  discussion,  a  variety  of 
approaches  to  the  identification  problem  exist.  For  example,  the  statistical 
parameter  estimation  procedures  assume  that  the  measured  modes  and  mass  and 
stiffness  matrices  all  contain  errors.  Some  procedures  assume  that  the 
measured  modes  are  exact  and  adjust  the  mass  and  stiffness  properties.  Others 
assume  Che  mass  matrix  to  be  exact,  orthogonalize  Che  measured  modes,  and  then 
adjust  the  stiffness  or  flexibility  matrix.  Apparently,  consensus  exists  only 
on  Che  need  to  adjust  analytical  stiffness  matrices. 


'  4 


-‘i 


1 ^Baruch,  M.  and  I.Y.  Bar  Itzhack,  "Optimal  Weighted  Orthogonalization  of 
Measured  Modes,"  AIAA  J. ,  Vol.  16,  No.  4,  Apr.  1978. 

^^Baruch,  M. ,  "Optimization  Procedure  Co  Correct  Stiffness  and  Flexibility 
Matrices  Using  Vibration  Tests,"  AIAA  J. ,  Vol.  16,  No.  11,  Nov.  1978. 

^^Berman,  A. ,  "Mass  Matrix  Correction  Using  an  Incomplete  Set  of  Measured 
Modes,"  AIAA  J. ,  Vol.  17,  No.  10,  Oct.  1979. 

^^Berman,  A.  and  E.J.  Nagy,  "Improvement  of  a  Large  Analytical  Model  Using 
Test  Data,"  AIAA  J. ,  Vol.  21,  No.  8,  Aug.  1983. 
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Analytical  dynamic  models  of  complex  structures  typically  involve  several 
hundred  degrees  of  freedom  (dof).  However,  at  best,  omy  a  few  dozen  measured 
modes  will  be  available.  Thus,  the  number  of  stiffness  coefficients  will 
greatly  exceed  the  number  of  equations  provided  by  the  measured  mode  data.  It 
is  reasonable  to  expect  that  a  more  accurate  identification  will  result  if  the 
ratio  of  stiffness  coefficients  to  available  equations  is  reduced.  This  can 
be  accomplished  by  supplementing  the  measured  mode  data  with  structural  connec¬ 
tivity  information.  Of  course,  by  doing  so,  we  assume  that  the  connectivity 
of  the  analytical  model  is  correct.  It  is  the  purpose  of  this  report  to 
describe  a  procedure  that  uses,  in  addition  to  mode  data,  structural  connec¬ 
tivity  information  to  optimally  adjust  deficient  stiffness  matrices. 

The  stiffness,  K,  matrix  adjustment  (KMA)  procedure  will  be  developed 
using  constrained  minimization  theory.  An  important  part  of  the  development 
is  the  selection  of  the  error  function  to  be  minimized.  In  Refs.  11  and  14, 
it  is  suggested  that  the  appropriate  error  function  should  involve  the  mass 
weighted  difference  between  the  analytical  and  adjusted  stiffness  matrix 
coefficients.  However,  in  complex  spacecraft  structures,  for  example,  non- 
structural  items  such  as  electronic  boxes,  batteries,  and  propellant  generally 
account  for  over  80  percent  of  the  total  mass.  Therefore,  it  is  inappropriate 
to  weigh  the  stiffness  coefficient  changes  by  any  function  involving  the  mass 
matrix.  In  addition,  the  error  function,  as  defined  in  Refs.  11  and  14,  can 
result  in  substantially  greater  percentage  changes  occurring  in  stiffness 
coefficients  with  numerically  small  values  than  in  coefficients  with  numer¬ 
ically  large  values.  This  bias  is  difficult  to  justify. 

Based  on  the  preceding  discussion,  it  must  be  concluded  that  the  error 
function  proposed  in  Refs.  11  and  14  is  not  suitable  for  our  purposes  and, 
therefore,  a  new  error  function  must  be  established.  It  is  proposed  that  the 
appropriate  error  function  be  independent  of  the  system  mass  properties  and 
stiffness  coefficient  magnitudes.  Furthermore,  minimization  of  the  error  func 
tion  should  minimize  the  percentage  change  to  each  stiffness  coefficient.  The 
optimally  adjusted  stiffness  matrix  can  then  be  obtained  by  minimizing  this 
error  function  suoject  to  symmetry  constraints,  connectivity  constraints,  and 
constraints  derived  from  the  system  eigenvalue  equation. 


2.  THEORETICAL  DEVELOPMENT 


The  preceding  discussion  indicated  that  it  is  reasonable  to  expect  a  more 
accurate  stiffness  matrix  adjustment  to  result  if,  in  addition  to  measured 
mode  data,  structural  connectivity  information  is  used.  This  can  be  accom¬ 
plished  by  insisting  that  all  coefficients  with  values  of  zero  in  the  original 
stiffness  matrix  also  have  values  of  zero  in  the  adjusted  stiffness  matrix. 
Mathematically,  this  can  be  achieved  by  insisting  that  the  adjusted  stiffness 
matrix  [K]  be  related  to  the  original  stiffness  matrix  [k]  by 

[K]  =  [k]0lY]  (1) 

where  [y]  is  a  matrix  to  be  determined,  and  the.  operator  0  defines  the 
following  element  operations: 


K.  . 
i-J 


k.  . 
ij 


(2) 


Therefore,  if  k. .  has  a  value  of  zero,  K. .  will  also  have  a  value  of  zero, 
ij  13 

Since  extensive  use  will  be  made  of  the  element-by-element  (scalar)  matrix 
multiplication  operator  0,  its  relevant  properties  are  presented  in  Appendix  A. 


To  minimize  unrealistic  changes  in  stiffness  coefficients,  the  error 
function  to  be  minimized  must  be  independent  of  the  stiffness  coefficient 
magnitudes  (see  Appendix  B).  This  can  be  accomplished  by  defining  the  error 
function  as 


where 


€  -  11(1]  -  (I|9[yI  1 1 


E  E  <i‘ij  -  i‘ij  nj: 

1=1  3=1 


I.  .  =  1 


=  0 


if 


if 


k.  .  ^  0 
ij 

k.  .  =  0 


(3) 


(4) 


The  adjusted  stiffness  matrix  can  now  be  obtained  by  estaolishing  the  [y] 
that  minimizes  G  and  satisfies  the  following  constraints: 


(5) 

(6) 

[M]  =  system  mass  matrix  (n,n)* 

[  4>]  =  system  mode  shapes  (n,m) 

2 

[oj^]  =  diagonal  matrix  of  circular  frequencies  squared  (m,m) 


-  [M][(t)][a)^]  +  ([k]0fY])[(t>]  =  [0] 
n 

[y]  -  [Y]"^  =  [0] 

where 


The  constraint  that  the  new  stiffness  matrix  be  consistent  with  the  mode 
shapes  and  natural  frequencies  is  introduced  by  Eq.  (5).  Inherent  in  this 
constraint  is  the  requirement  that  the  mode  shapes  exhibit  the  orthogonality 
associated  with  normal  modes,  i.e. 

[<{']'^[Ml[<{)]  =  [I]  (7) 

The  requirement  that  the  new  stiffness  matrix  be  symmetric  is  introduced  by 
Eq.  (6) . 


The  method  of  Lagrange  Multipliers^^  will  be  used  to  incorporate  in  the 
minimization  of  e  the  constraints  defined  by  Eqs.  (5)  and  (6).  We  begin  by 


forming  the  Lagrange 

Function 

L 

n 

m 

/  n  \ 

L  = 

E 

i=l 

E  • 

j-i 

(  hi)  -  ^ 

n 

n 

>  E 

i=l 

'Mj  -  i'ji> 

(8) 

"'Designates  dimension  of  matrix 

'■^Hadley,  G. ,  Nonlinear  and  Dynamic  Programming,  Addison-Wes ley  Publishing 
Company,  Inc.,  Reading,  Mass.,  1964. 


Jiiere  A.  .  and  u.  .  are  the  La^ran^e  Multipliers.  The  quantitv  E  represents 
ijij- 

the  first  term  of  Eq.  (5),  which  is  not  a  function  of  and,  therefore, 

does  not  need  to  be  explicitly  defined.  Proceeding,  we  take  the  partial 
derivative  of  L  with  respect  to  each  Y--*  These  derivatives  are  set  equal 

'ij 

to  zero  to  obtain  a  set  of  nxn  equations  that  the  must  satisfy  for  L 
to  be  a  minimum.  Expressing  these  equations  in  matrix  notation,  we  obtain 


-2([I]  -  [y])  +  (k]@([X]((t)]^)  +  [y]  =  [0]  (9) 

Equations  (5),  (6),  and  (9)  represent  the  complete  system  of  equations 

T 

needed  to  determine  the  unknowns  (X],  [ y] ,  and  [y].  Since  [y]  =  -[y]  , 

[y]  can  be  eliminated  by  adding  Eq.  (9)  and  its  transpose 

-4([I]  -  fy])  +  [k]@([X][(t)]'^  +  I(|)][X]^)  =  [0]  (10) 

Next,  we  pre-e lement-by-eleraent  multiply  Eq.  (10)  oy  l/4[k]  and  rearrange 
terms  to  obtain 

M@[y]  =  Ik]  -  1/4  14>]©([X]U]^  +  (<1)][X]'^)  (11) 

where 

[<t]  =  (k]@[k]  (12) 

Substituting  Eq.  (11)  into  Eq.  (5),  we  obtain 

[A]  +  {[$]@([X]f(t)]'^)}(4)]  +  {[<t]0([4)](Xj'^)}[(i)]  =  [0]  (13) 

where 

[A]  =  4([M]l(|)][w^]  -  lk]l()>])  (14) 

n 

Equation  (13)  can  now  be  used  .to  establish  [  X]  which,  when  substituted  into 
Eq.  (11),  yields  the  desired  stiffness  matrix  [see  Eq.  (1)1. 

3y  performing  the  operations  indicated  in  Eq.  (13),  and  taking  advantage 
of  relationship  5  presented  in  Appendix  A,  it  can  be  shown  that 


syminecric  correction  matrix  orthogonal ization  procedure  described  in  Ref.  16  will 
be  used,  i.e. 

U*"]  =  (23) 

where  [  (J*™]  are  the  measured  modes  and  [4>'’l  are  the  orthogonalized  modes. 

Only  the  first  five  modes  were  included  in  the  orthogonalization,  whereas  all 
eight  modes  contributed  in  the  response  calculations.  This  simulates  actual 
test  conditions  where  modes  that  are  not  measured  contaminate  the  measured 
modes.  The  two  orthogonalized  modes  used  in  the  subsequent  stiffness  matrix 
adjustments  are  compared  in  Table  9  to  the  "test"  modes  and  the  true  normal 
modes. 

The  adjusted  stiffness  matrix  obtained  using  the  first  orthogonalized 
"test"  mode  is  presented  in  Table  10.  Comparing  these  coefficients  to  the 
exact  stiffness  matrix  (Table  1),  it  can  be  observed  that  many  of  the 
"analytical"  stiffness  matrix  coefficients  have  improved  substantially, 
particularly  the  diagonal  terms.  By  comparing  the  adjusted  matrix  coef¬ 
ficients  to  those  obtained  with  a  single  normal  mode  (Table  4),  it  can  be 
observed  that  they  are  comparable.  In  addition,  connectivity  has  been 
preserved. 

The  adjusted  stiffness  matrix  obtained  with  two  orthogonalized  "test" 
modes  is  presented  in  Table  11.  By  comparing  the  adjusted  stiffness  matrix 
coefficients  to  those  in  Table  5,  we  observe  comparable  improvement  to  that 
obtained  with  two  normal  modes.  In  addition,  as  has  been  the  case  for  all 
test  problems,  connectivity  has  been  preserved.  Thus,  with  a  limited  number 
of  mode  shapes,  knowledge  of  the  structural  connectivity,  and  a  physically 
realistic  minimization  function,  it  is  possible  to  improve  a  deficient  stiff¬ 
ness  matrix. 


^^Targoff,  W.  P. ,  "Orthogonality  Check  and  Corre  tion  of  Measured  Modes," 
AIAA  J. ,  Vol.  14,  No.  2,  Feb.  1976. 
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Table  7.  Comparison  of  Natural  Frequencies 


Natural  Frequencies  (rad/sec) 

Corrupted  Stiffness 

Adjusted 

Stiffness 

Matrix 

Exact  Stiffness 

Matrix 

1  Mode 

2  Modes 

3  Modes 

Matrix 

24.769 

30.665* 

30.665* 

30.665* 

30.665 

38.776 

31.737 

31.711* 

31.711* 

31.711 

38.859 

31.850 

31.763 

31.763* 

31.763 

31.529 

33.694 

32.362 

32.362 

32.362 

41.783 

37.317 

34.194 

34.193 

34.193 

42.015 

39.443 

35.561 

35.560 

35.560 

44.811 

41.393 

38.789 

38.789 

38.789 

54.S41 

50.590 

46.412 

41.888 

41.888 

*  Modes  used  to  adjust  corrupted  stiffness  matrix 


Table  6.  Adjusted -Stiffness  Matrix  Coefficients — Three  Normal  Modes 


2 

3 

4 

5 

6 

7 

8 

■ 

1.5 

-1.5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1011.5 

-10.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1110.0 

0.0 

-100.0 

0.0 

0.0 

0.0 

1100.0 

-100.0 

-100.0 

0.0 

0.0 

SYM 

1100.0 

0.0 

0.0 

0.0 

1112.0 

-10.0 

-2.0 

1011.5 

-1.5 

8 

3.5 
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The  stiffness  matrix  coefficients  identified  with  three  normal  modes  are 
presented  in  Table  6.  By  comparing  these  coefficients  to  the  exact  values 
presented  in  Table  1,  it  can  be  observed  that  they  are  identical.  Thus,  with 
a  subset  of  the  system  normal  modes,  and  knowledge  of  the  structure's  con¬ 
nectivity,  the  exact  stiffness  matrix  coefficients  have  been  identified. 

The  natural  frequencies  associated  with  each  of  the  adjusted  stiffness 
matrices  are  compared  in  Table  7  to  the  exact  values.  The  correspondence 
between  mode  shapes,  and  associated  frequencies,  was  established  by  performing 
a  cross-orthogonality  check  between  the  true  normal  modes  and  the  modes  asso¬ 
ciated  with  each  of  the  adjusted  stiffness  matrices.  As  can  be  ascertained 
from  the  table,  the  natural  frequencies  of  the  modes  used  in  the  identifica¬ 
tion  are  reproduced  exactly  by  the  adjusted  model.  A  comparison  of  corres¬ 
ponding  raoii», shapes  also  showed  exact  agreement.  Furthermore,  for  the  one- 
and  two-mode  cases,  the  modes  and  frequencies  that  were  not  used  in  the  identi¬ 
fication  also  are  in  closer  agreement  with  the  true  values.  The  three-mode 
case,  as  expected,  provides  exact  agreement  in  all  modes. 

To  be  of  practical  use,  the  KMA  procedure  must  yield  reasonable  results 
when  imperfect  test-derived  modes  are  used.  To  study  the  procedure's  sen¬ 
sitivity  to  imperfect  data,  the  applicable  test  cases  were  repeated  using 
simulated  test  modes.  The  test  modes  were  obtained  by* establishing  the  closed 
form  quadrature  response  of  the  structure  to  multiple  sinusoidally  varying 
forces  (A  critical  damping  ratio  of  0.01  was  assigned  each  mode.).  The 
multiple  force  levels  were  adjusted  to  obtained  "measured"  modes  whose  mutual 
contamination  did  not  exceed  approximately  10  percent.  The  unit-normalized 
generalized  mass  matrix  for  the  first  five  modes  is  presented  in  Table  8.  The 
level  of  contamination  is  what  can  be  expected  from  a  properly  performed  mode 
survey  test  of  a  complex  structure. 

Before  proceeding  with  the  identification,  it  is  necessary  for  the  "test" 
modes  to  be  analytically  orthogonalized.  For  the  purposes  of  this  study,  the 
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Table  4.  Adjusted  Stiffness  Matrix  Coefficients — One  Normal  Mode 


2 

3 

4 

5 

6 

7 

8 

■ 

1.7 

-2.1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1013.5 

-10.1 

0.0 

0.0 

0.0 

0.0 

0.0 

1275.8 

0.0 

-198.6 

0.0 

0.0 

0.0 

1237.4 

-178.6 

-198.5 

0.0 

0.0 

1237.5 

0.0 

0.0 

0.0 

1279.7 

-10.1 

-4.1 

1016.1 

-2.0 

8 

5.1 

Table  5.  Adjusted  Stiffness  Matrix  Coefficients — Two  Normal  Modes 


Initially,  it  appears  that  all  16  stiffness  coefficients  can  be  identified 
from  two  normal  modes.  However,  as  illustrated  in  Fig.  1,  dof  6,  7,  and  8 
involve  a  load  path  indeterminacy  of  three.  Thus,  to  obtain  exact  identi¬ 
fication,  three  normal  modes  are  needed.  If  less  than  three  normal  modes  are 
used,  an  adjusted  stiffness  matrix  will  be  obtained  that  minimizes  the  error 
function  [Eq.  (3)]  and  satisfies  the  symmetry  [Eq.  (6)),  eigenproblem 
[Eq.  (5)],  and  connectivity  [Eq.  (1)]  constraints.  To  demonstrate  this,  the 
"analytical"  stiffness  matrix  was  first  adjusted  using  only  one,  then  two,  and 
finally  three  normal  modes. 

The  adjusted  stiffness  matrix  obtained  using  one  normal  mode  is  presented 
in  Table  4.  Comparing  this  matrix  to  the  exact  stiffness  matrix,  shown  in 
Table  1,  we  observe  that  many  of  the  "analytical"  stiffness  matrix  coefficients 
have  significantly  approached  the  true  values.  In  addition,  connectivity  has 
been  preserved. 

The  adjusted  stiffness  matrix  obtained  using  the  first  two  normal  modes  of 
the  structure  is  presented  in  Table  5.  Again,  connectivity  has  been  preserved. 
In  addition,  the  load  path  stiffnesses  associated  with  dof  1  through  5  are  in 
exact  agreement  with  the  coefficients  presented  in  Table  1,  Furthermore,  the 
terras  associated  with  dof  6  and  7  have  also  improved  substantially. 

This  two-mode  case  illustrates  an  interesting  feature  of  the  KMA  procedure. 
The  procedure  will  identify  exactly  those  stiffness  coefficients  associated 
with  the  part  of  the  structure  for  which  the  load  path  indeterminacy  does  not 
exceed  the  number  of  modes  used.  The  stiffness  coefficients  associated  with 
the  part  of  the  structure  whose  redundancy  exceeds  the  number  of  modes  used 
will  be  adjusted  to  minimize  the  error  function  of  Eq.  (3)  and  satisfy  the 
connectivity,  eigenproblem,  and  symmetry  constraints  [Eqs.  (1),  (5),  and  (6)]. 

As  discussed  previously,  three  normal  modes  would  provide  the  needed 
conditions  to  identify  the  stiffness  matrix  coefficients  exactly.  Note  that 
regardless  of  whether  three,  four,  or  all  eight  modes  are  used,  [a]  +  [S] 
will  yield  only  16  nonzero  eigenvalues,  and  the  resulting,  stiffness  matrices 
will  all  be  identical. 
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2 


3 


5 


6 


8 


I 


2.0  -2.0  0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1512.0  -10.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1710.0 

0.0 

-200.0 

0.0 

0.0 

0.0 

850.0 

-200.0 

-200.0 

0.0 

0.0 

850.0 

0.0 

0.0 

0.0 

SYM 

1714.0 

-10.0 

-4.0 

1512.0 

-2.0 

6.0 

3.  DEMONSTRATION  OF  PROCEDURE 


The  KMA  procedure  will  be  demonstrated  by  numerical  simulation  of  a  test 
problem.  The  procedure  will  be  used  to  adjust  the  corrupted  stiffness  matrix 
of  an  8  dof  analytical  structure.  The  adjustments  will  first  be  performed 
using  the  normal  modes  of  the  system.  The  problem  will  then  be  repeated  using 
simulated  test  modes. 

A  schematic  representation  of  the  structure  is  shown  in  Fig.  1,  where  the 
squares  represent  dof  and  the  springs  represent  load  paths.  The  stiffness 
matrix  of  the  structure  is  presented  in  Table  1,  and  the  diagonal  terms  of  the 
diagonal  mass  matrix  are  shown  in  Table  2.  This  structure  represents  a  severe 
test  case  because  of  the  large  relative  differences  in  the  magnitudes  of  some 
of  the  stiffness  matrix  coefficients.  In  addition,  all  eight  natural  frequen¬ 
cies  are  within  27  percent  of  each  other,  and  the  first  four  natural  frequen¬ 
cies  are  within  5  percent  of  each  other. 

The  stiffness  matrix  that  represents  the  best  "analytical"  model  of  the 
structure  is  presented  in  Table  3.  This  matrix  was  obtained  by  modifying  the 
stiffness  matrix  coefficients  of  Table  1.  Some  values  were  increased  by  a 
maximum  of  100  percent,  others  were  not  changed,  and  two  diagonal  terms  were 
decreased  by  nearly  30  percent.  Comparing  the  two  matrices,  it  can  be 
observed  that  the  "analytical"  stiffness  matrix  does  not  resemble  the  true 
stiffness  matrix,  except  for  connectivity. 

Before  proceeding  with  the  matrix  adjustment,  the  number  of  independent 
coefficients  in  the  "analytical"  stiffness  matrix  will  be  established.  As 
discussed  previously,  this  can  be  accomplished  by  adding  the  number  of 
diagonal  terms  in  the  matrix  and  the  number  of  nonzero  coefficients  on  one 
side  of  the  diagonal.  For  the  test  problem,  16  independent  terms  are  obtained 


[4)]^[a]  +  [6])[iH{p} 


(21) 


where  [(2]  is  a  diagonal  matrix  whose  diagonal  terms  are  the  nonzero  eigen¬ 
values  of  [ot]  +  [6].  Using  Eq.  (21)  to  solve  for  {p}  and  substituting  into 
Eq.  (20),  we  finally  obtain 

{X}  =  (22) 

We  can  now  construct  [X]  [see  Eq.  (17)]  and  establish  [K]  using  Eqs.  (11) 
and  (1). 

The  transformation  defined  by  Eq.  (20)  yields  a  solution  [Eq.  (22)]  that 
is  applicable  when  the  number  of  constraints  defined  by  Eq.  (5)  exceeds  the 
number  of  independent  stiffness  coefficients  available  for  adjustment.  If 
normal  modes  are  used,  the  constraints  will  be  consistent  with  each  other  and 
the  true  properties  of  the  structure.  Thus,  the  inclusion  of  additional  modes 
past  a  certain  threshold,  will  not  alter  the  identified  stiffness  coefficients 

This  will  not  be  the  case  if  imperfect  test  modes  are  used.  Because  of 
measurement  error,  empirical  modes  will  not  necessarily  be  perfectly  con¬ 
sistent  with  each  other  and  the  true  properties  of  the  structure.  However, 

Eq.  (22)  will  still  provide  a  solution.  Discussion  of  this  feature  is  beyond 
the  intended  scope  of  this  presentation.  For  the  purposes  of  the  present  dis¬ 
cussion,  if  we  are  dealing  with  test  modes,  we  shall  restrict  our  attention  to 
the  classical  constrained  minimization  problem  in  which  the  number  of  con¬ 
straints  does  not  exceed  the  number  of  coefficients  available  for  adjustment. 


where 


[H]..  =  - 


.th 

J 


column  of 


{D^}^  *  i*"^  column  of  (<t]^[4'^] 


Equation  (15)  can  now  be  used  to  establish  {X}.  Note  that  the  rank  of 
[a]  +  [6]  will  not  exceed  a  number  equal  to  the  number  of  diagonal  terms  in  [k] 
plus  all  the  nonzero  coefficients  on  one  side  of  the  diagonal.  The  exact  rank 
of  [a]  +  [S]  will  depend  on  the  structural  connectivity  and  the  number  of  inde¬ 
pendent  constraints  defined  by  Eq.  (5).  For  many  problems  of  practical  interest 
(a]  +  [6]  will  be  singular.  In  addition,  [ot]  and  [3]  are  symmetric  and 
[a]  +  (3]  is  indefinite;  i.e.,  the  eigenvalues  of  [a]  +  [31  can  be  negative, 
zero,  and/or  positive. 


Equation  (15)  can  be  solved  by  first  defining  the  following  transformation: 

{X}  =  fil^Kp}  (20) 

where  the  columns  of  [ip]  are  the  eigenvectors  associated  with  the  nonzero  eigen¬ 
values  of  [a]  +  [3].  Substituting  Eq.  (20)  into  Eq.  (15)  and  premultiplying  by 
T 

[ip]  ,  we  obtain 


{a}  =  ([a]  -K  [61){x} 

where  the  elements  of  [a]  and  [X]  have  been  written  as  column  vectors  {a} 
and  {X},  respectively.  For  example 


11 

Ai2 

21 

A22 

31 

A32 

{a}  = 


[X]  = 


11 

\2 

21 

^22 

31 

^2 

with  {a}  and  {X}  defined  by  Eqs.  (16)  and  (17),  it  can  be  shoim  that  [a] 
and  [B]  are  as  follows: 


[G^]  I 


i  T  ^  i  ^  i 

where  (G  ]  =  -[(t)]  [$  ][<{)]  and  J  is  a  diagonal  matrix  whose  diagonal  terms 


are  the  i  row  of  [<I>]  [see  Eq.  (12)1,  and 


Table  9.  Comparison  of  Mode  Shapes 


"Test"  Mode  Shapes 

OrChogonalized  "Test" 

Normal  Mode 

Mode  Shapes 

Shapes 

0.029 

0.112 

0.144 

0.015 

0.042 

0.054 

0.374 

0.356 

0.360 

0.595 

0.609 

0.606 

{(|)}l  0.604 

0.600 

0.605 

0.356 

0.365 

0.362 

0.117 

0.076 

0.062 

0.551 

0.522 

0.505 

2.753 

2.774 

2.776 

0.909 

0.917 

0.915 

0.098 

0.138 

0.127 

-0.080 

-0.078 

-0.088 

{(J)}2  0.927 

0.050 

0.041 

-0.127 

-0.123 

-0.124 

-0.365 

-0.330 

-0.337 

-0.539 

-0.498 

-0.506 
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Table  10.  Adjusted  Stiffness  Matrix  Coefficients — One  "Test"  Mode 


4 .  SUMMARY 


A  procedure  has  been  introduced  that  uses,  in  addition  to  mode  data, 
structural  connectivity  information  to  optimally  adjust  deficient  stiffness 
matrices.  The  procedure  was  developed  using  constrained  minimization  theory. 
The  minimization  error  function  was  formulated  such  that  the  resulting  changes 
to  stiffness  coefficients  are  a  minimum.  The  resulting  procedure  retains 
the  physical  configuration  of  the  analytical  model,  and  the  adjusted  model 
exactly  reproduces  the  modes  used  in  the  identification. 

The  stiffness  matrix  adjustment  (KMA)  procedure  was  demonstrated  by 
numerical  simulation  of  a  test  problem.  The  procedure  was  used  to  adjust  the 
corrupted  stiffness  matrix  of  an  8  dof  analytical  structure.  The  adjustments 
’were  first  performed  using  the  system  normal  modes,  of  which  only  three  were 
needed  for  exact  identification.  The  test  problem  was  then  repeated,  using 
simulated  test  modes,  with  excellent  results. 
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APPENDIX  A 


PROPERTIES  OF  ELEMENT-BY-ELEMENT 
MATRIX  MULTIPLICATION  OPERATOR 

The  derivation  of  the  KMA  procedure  used  an  element-by-element  matrix 
multiplication  operator,  0.  Relevant  properties  of  the  operator  are  presented 
below: 


[c]  =  (a]@[b] 


defines  c . .  =  a . .b . . 
ij  iJ 


[c]  =  [d]([a]0[b)) 


II 

implUs  c.  .  - 


(a]0[b]  =  (b)0[a] 


{[d]([a]0(b])}^  =  ((a]^©lb]’^)[d]^ 


(a]0[b]  =  ^  [a]  [b] 

nxn  nxn  j=l 


•ik  ■  *jk  ■  i  '■te  ■  "jk 


=  0  for  i'  ^  i 


Ok--k’  k^'  -v  -  •-  JV  J  m  — *i 


=  0  for  i  k 


APPENDIX  B 


ADJUSTMENT  PROCEDURE  WITH  STIFFNESS 
COEFFICIENT  MAGNITUDE  WEIGHING 

la  the  theoretical  development  of  the  KMA  procedure,  it  was  asserted 
that,  to  minimize  unrealistic  changes  in  stiffness  coefficients,  the  error 
function  must  be  independent  of  the  stiffness  coefficient  magnitudes.  This 
was  accomplished  by  defining  the  error  function  c  as 


e  =  llfl]  -  [I]  0  [y]I‘ 


(B-1) 


where  [l]  was  obtained  by  replacing  the  nonzero  coefficients  of  [k]  with  1.0s. 

It  is  a  relatively  simple  task  to  explore  the  impact  of  incorporating, 
in  the  error  function,  the  stiffness  coefficient  magnitudes.  This  can  be 
accomplished  by  redefining  the  error  function  as 

e  =■  ll{k]  -  [kl  0  [y)ii 


(B-2) 


The  symmetry  and  eigenproblem  constraints  are  introduced  in  the  minimization 
of  the  above  error  function,  as  before,  using  Lagrange  Multipliers,  to  obtain 


-2([k]0[k]  -  [k]0(k]0(Y])  +  (k)0((Xll<))]")  +  [y]  *  [0] 


(B-3) 


We  eliminate  [y]  by  adding  Eq.  (B-3)  and  its  transpose 


-4([k]@(k]  -  [k]0[k]0[Y])  +  [k]0((X][(|)]^+  (♦IIX)’^)  -  [0] 


(B-4) 


Proceeding,  we  define  a  matrix  [k]  such  that 


k.  .  =■  1/k.  .  if  k.  .  0 

ij  iJ 


(B-5) 


if  k.  .  =■  0 
ij 


— I  fi  II  iTttS- iTrf -rt i"i  i^t-  **  ft , 
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Mext,  we  pre-eleraent-by-element  multiply  Eq.  (B-4)  Sy  1/4  [k]  and  rearrange 
terms  to  obtain 

[k]@[Y]  =  [k]  -  1/4  [I]@([X][<{)]'^  +  [(t>][X]'^)  (3-6) 

where  all  terms  are  as  defined  in  Section  2. 

Comparing  Eq.  (B-6)  to  Eq.  (11),  we  note  that  they  are  identical  except 

A 

that  [l]  has  replaced  [$].  Therefore,  the  solution  as  defined  by  Eq.  (22) 

A 

is  applicable,  and  the  only  needed  change  is  the  substitution  of  [I]  for  [ $] . 

The  analytical  test  structure  corrupted  stiffness  matrix  (Table  3)  was 
adjusted  using  one  and  then  two  normal  modes.  As  with  the  KMA  procedure, 
three  modes  provide  exact  identification,  and  both  formulations  yield  exact 
agreement  for  coefficients  defined  by  constraints  only.  The  adjusted  matrix 
coefficient,  presented  in  Table  B-1  (one  mode)  and  Table  B-2  (two  modes)  can 
be  compared  to  the  corresponding  values  obtained  with  the  KMA  procedure 
(Tables  4  and  5)  and  to  the  exact  values  (Table  1).  As  can  be  ascertained, 

Che  coefficients  identified  with  Che  KMA  procedure  are  in  considerably  better 
agreement  with  Che  exact  values  Chan  Che  coefficients  obtained  using  Che  error 
function  defined  by  Eq.  (B-2). 

Particularly  troublesome,  with  the  error  function  defined  by  Eq.  (B-2), 
are  Che  large  changes  that  occur  in  stiffness  coefficients  with  small  relative 
numerical  values.  For  example,  the  diagonal  term  associated  with  dof  8  was 
modified  from  a  value  of  6.0  (Table  3)  to  a  value  of  116  [Table  (B-2)].  This 
should  be  compared  to  the  adjustment  made  by  the  KMA  procedure  which  changed 
the  coefficient  from  the  value  of  6.0  to  a  value  of  4.3.  The  true  value  for 
this  particular  coefficient  is  3.5  (Table  1). 

The  poor  performance  of  the  error  function  defined  by  Eq.  (B-2)  is  due 
to  the  inclusion  of  the  stiffness  coefficient  magnitudes.  The  minimization  of 
an  error  function,  such  as  defined  by  Eq.  (B-2),  will  result  in  small  per¬ 
centage  changes  to  large  magnitude  coefficients  and  relatively  large  percent¬ 
age  changes  to  small  magnitude  coefficients.  Therefore,  to  minimize  unreal¬ 
istic  changes  in  stiffness  coefficients,  the  error  function  should  be  inde¬ 
pendent  of  the  stiffness  coefficient  magnitudes. 
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Table  B-1.  Adjusted  Stiffness  Matrix  Coefficients — One  Normal  Mode 


Table  B-2.  Adjusted  Stiffness  Matrix  Coefficients — Two  Normal  Modes 
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1.5 

-1.5  0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9 

1011.5  -10.0 

0.0 

0.0 

0.0 

0.0 

0.0 

9 

1110.0 

0.0 

-100.0 

0.0 

0.0 

0.0 

9 

1100.0 

-100.0 

-100.0 

0.0 

0.0 

5 

SYM 

1100.0 

0.0 

0.0 

0.0 

6 

1278.0 

135.0 

-139.0 

7 

1137.0 

-121.0 

8 

116.0 

NOMENCLATURE 
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ij 

(Kl 

K.  . 
ij 
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see  Eq.  (14) 

vector  of  [A]  elements  [see  Eq.  (16)] 
defined  by  Eq.  (18) 
defined  by  Eq.  (19) 

i^*^  column  of  [<^]^  1$'^] 

matrix  of  stiffness  matrix  adjustment  coefficients 
element  ij  of  [yl 

error  function 

submatrix  i  of  [a] 
submatrix  ij  of  [g] 

identity  matrix 

see  Eq.  (4) 

analytical  stiffness  matrix 
element  ij  of  [k] 

adjusted  analytical  stiffness  matrix 
element  ij  of  [k] 

Lagrange  function 
Lagrange  Multiplier 


1 


•'1 


.■J 


1 
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NOMENCLATURE  (Continued) 


matrix  of  Lagrange  Multipliers 

vector  of  [X]  elements  (see  Eq.  (17)] 
mass  matrix 
Lagrange  Multiplier 

Matrix  of  Lagrange  Multipliers 

element-by-element  matrix  multiplication  operator 

transformed  {X}  [see  Eq.  (20)] 
matrix  of  normal  mode  vectors 

i^^  normal  mode  vector 

column  of  [<|)]^ 
element  ij  of  (4)] 

matrix  of  measured  mode  vectors 

matrix  of  analytically  orthogonalized  measured  mode  vectors 
[k]  0  [k] 

diagonal  matrix  whose  diagonal  elements  are  the  i^*^  row  of  ($] 
eigenvectors  of  [a]  +  (3)  associated  with  [il] 

diagonal  matrix  of  circular  frequencies  squared 

nonzero  eigenvalues  of  [a]  *  18] 


null  matrix 


